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Steen Hannestad 
NORDITA, Blegdamsvej 17, 2100 Copenhagen, Denmark 



H. -Thomas Janka 

Max- Planck- Institut fur Astrophysik, Karl-Schwarzs child- Str. 1, 85740 Garching, Germany 

Georg G. Raffelt 

Max- Planck- Institut fur Physik (Werner-Heisenberg-Institut) , Fohringer Ring 6, 80805 Miinchen, Germany 
and TECHNION, Israel Institute of Technology, Haifa 32000, Israel 



o 
o 
o 

(N 

< 

(N 
> 
(N 
^J- 
(N 
(N 

On 
On 



9 L,: 
6 

CO 



X 



Giinter Sigl 

Departement d'Astrophysique Relativiste et de Cosmologie, 
UMR 8629 du CNRS, Observatoire de Paris, 5 place J. Janssen, 92195 Meudon Cedex, France 

(8 December 1999, revised 19 May 2000) 

We study if the neutrino mixing parameters suggested by the atmospheric neutrino anomaly imply 
chemical equilibrium between fi- and r-flavored leptons in a supernova (SN) core. The initial flavor- 
conversion rate would indeed be fast if the v^-Vr- mixing angle were not suppressed by second-order 
refractive effects. The neutrino diffusion coefficients are different for z/ M , v T and v T so that 
neutrino transport will create a net \x and r lepton number density. This will typically lead to 
a situation where the usual first-order refractive effects dominate, further suppressing the rate of 
flavor conversion. Altogether, neutrino refraction has the nontrivial consequence of guaranteeing 
the separate conservation of e, fi, and r lepton number in a SN core on the infall and cooling time 
scales, even when neutrino mixing angles are large. 



PACS numbers: 97.60.Bw, 14.60.Pq 



I. INTRODUCTION 

The evidence for neutrino flavor oscillations from at- 
mospheric 01 and solar neutrinos |Q] as well as the LSND 
experiment |3j is now so compelling that the debate in 
neutrino physics has fundamentally changed. The cur- 
rent experimental results point to a few very specific re- 
gions in the parameter space of mass differences and mix- 
ing angles. Therefore, the experimental task at hand is 
to verify the proposed solutions by independent means 
such as long-baseline oscillation experiments. From the 
astrophysical perspective, it is no longer an academic ex- 
ercise to investigate systematically if neutrino mixings 
with these specific parameters lead to observable effects 
other than explaining the solar neutrino deficit. 

One consequence of the apparent violation of e, fi, 
and t flavor lepton number is that the leptons of a per- 
fectly thermalized system are no longer characterized by 
six independent chemical potentials, but only by one 
for the charged leptons and one for the neutrinos, i.e. 
A*e = Mm = Mr and fi v<! = n v = jj, Ur . In practice one 
knows of only two types of environment where neutrinos 
achieve thermal equilibrium and where this effect could 
be important, the hot early universe and the inner cores 
of collapsed stars. In the standard picture of the early 
universe, the neutrino chemical potentials are very small 
so that flavor equilibration would not make much of a 
difference except for the case of oscillations between ac- 



tive and sterile neutrinos 0, a possibility that we will 
mostly ignore. 

In a supernova (SN) core, on the other hand, the elec- 
tron lepton number of the progenitor star's iron core is 
trapped, leading to highly degenerate e and v e distribu- 
tions. If the trapped electron lepton number were shared 
among all flavors, the impact on the effective equation 
of state would be large. Equilibration on the infall time 
scale of a few hundred milliseconds could well be disas- 
trous for the SN explosion mechanism. However, in a 
previous paper || two of us found that even for max- 
imal neutrino mixing, equilibrium between v e and the 
other flavors is not achieved on the time scale of a few 
seconds unless 5m 2 exceeds about 10 5 eV 2 . The rate of 
flavor conversion is very slow because the mixing angle is 
strongly suppressed by neutrino refractive effects. Since 
the current indications for neutrino oscillations point to 
much smaller mass differences, the electron lepton num- 
ber in a SN core is almost perfectly conserved except for 
diffusive or convective transfer to the stellar surface. 

We presently study the corresponding conversion be- 
tween and u T , which has not yet been investigated. It 
is usually assumed that the chemical potentials for these 
flavors vanish in a SN core so that their distribution func- 
tions are equal, i.e. that flavor equilibrium exists from the 
start. However, the muon mass of 106 MeV is not very 
large, considering that typical temperatures can exceed 
T = 30 MeV and that the average thermal energy of a 
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relativistic particle is around 3T. Beta equilibrium by re- 
actions of the form i/^ + n <-> p + implies the condition 
A/i = ii n — \Lp = fJi^ — fJ-Vai which is familiar from the elec- 
tron flavor. The exact value of A/x depends sensitively 
on the equation of state; a typical range is 50-100 MeV. 
Since initially the trapped p, lepton number is zero, and 
taking A/i = 50 MeV and T — 30 MeV as an example, 
we find /i^ w 32 MeV and fi v w — 18 MeV so that there 
is a significant excess of pT over p + , compensated by an 
equal excess of ZA, over i/„. On the other hand, the large 
value m r = 1777 MeV completely suppresses the pres- 
ence of r leptons. Together with the absence of trapped 
t lepton number this implies p Vr = 0. The initial ther- 
mal distributions are quite different between the p and r 
flavors! 

Apart from these initial differences, it was recently 
recognized that the opacities for v^, v^, v T and v T are 
all different from each other if nucleon recoil effects and 
muonic beta reactions are taken into account, leading to 
a transient build-up of p and t lepton number as neu- 
trinos diffuse out of the star ||. Again, it is important 
to understand if the p and r lepton numbers are locally 
conserved on the diffusion time scale or if the condition 
fi v = \l Vt is enforced by flavor-conversion processes. 

The SuperKamiokande measurement of the atmo- 
spheric neutrino anomaly M suggests maximal v^-Vr- 
mixing and Sm 2 roughly between 10 -3 and 10~ 2 eV 2 . 
We now study if these mixing parameters lead to chemi- 
cal equilibrium between and v T on the diffusion time 
scale of a few seconds. 



II. FLAVOR CONVERSION 

Chemical equilibrium between charged leptons £ and 
neutrinos vi of a given flavor is established by beta pro- 
cesses of the form t~ + p <-> n + V£. They are "infinitely 
fast" relative to all other time scales of interest. Chemi- 
cal equilibrium between p- and r-flavored leptons, on the 
other hand, is achieved by neutrino oscillations and by 
collisions which break the coherence of mixed neutrino 
states. The evolution of the f M -i/ r -system is best de- 
scribed in terms of the 2x2 density matrices p p for each 
momentum p. In the flavor basis, the diagonal entries 
of p p are the usual and v T occupation numbers, re- 
spectively, while the off-diagonal entries represent phases 
produced by oscillations. The evolution of p p , and of p p 
for anti-neutrinos, is governed by a generalized Boltz- 
mann collision equation which simultaneously includes 
the effects of oscillations and collisions |7],|8| . 

For our simple estimates, however, it will suffice to use 
a "pedestrian" method where the neutrino ensemble is 
represented by a single mode with momentum p or en- 
ergy E « |p|. Without collisions, the flavor content of 
this mode would oscillate forever. Interactions with the 
medium, however, destroy the coherence between the fla- 
vor components, leading eventually to an equal, incoher- 



ent mixture. Maximally mixed neutrinos approach this 
state of flavor equilibrium with a rate |9(| 



2 D 
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where 
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= 1.7 x 10~ n eV AsE^ 1 
= 2.5 x 10 4 s" 1 A 3 E^ 1 
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is the oscillation frequency in vacuum. Here, A3 = 
(5m 2 /10 -3 eV 2 represents the lower end of the mass range 
implied by SuperKamiokande and E 30 = E/30 MeV 
where 30 MeV is a typical temperature in a SN core. 
Further, SV is the energy difference between and v T 
of equal momenta caused by the medium, i.e. SV is the 
difference of the medium's weak potential for our two 
neutrino flavors. 

Finally, D is the damping or decoherence rate, i.e. the 
rate by which interactions with the medium "measure" 
the flavor content of a mixed neutrino state. Typically D 
is of the order of the neutrino collision rate, but the exact 
relationship between the two quantities is not trivial. For 
example, in a situation where one of the neutrino flavors 
scatters with a rate r co n while the other is sterile, one 
finds D = r co n/2 ||. Equation ([!]) applies in the "strong 
damping limit" defined by D 3> w V ac, a condition which 
is satisfied in our scenario. 

Equation (Q) is easily interpreted in two limiting cases. 
For SV = 0, the flavor content of a given state oscillates 
as i[l + cos(a; vac t)]. The oscillations are interrupted by 
those collisions which "measure" the difference between 



Vu and v T 



If this "measurement rate" or "decoher- 



ence rate" D = t^ 1 is much larger than w vac , the flavor 
oscillation is interrupted when w vac t C 1 so that the fla- 
vor content evolves as (w vac i/2) 2 until it is interrupted. 
This happens with a rate r^ 1 so that the rate of flavor 
conversion must scale as (w vac T]j) 2 t^ 1 or 
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D 
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In this case the flavor conversion rate decreases with in- 
creasing D, and vanishes for infinite D. This situation 
is known as the "Quantum Zeno Paradox" or "Watched 
Pot Effect" ||: The neutrino remains "frozen" in its fla- 
vor state because it is frequently "watched" or measured 
to be in this state by the interactions with the medium. 

A more familiar limiting case obtains when \SV\ D 
so that D 2 in the denominator of Eq. (|l|) can be ne- 
glected. Since D 3> uj viLC by assumption, we also have 
\SV\ 3> w ra c, implying that the oscillation frequency is 
SV instead of w vac because the energy difference between 

and v T of equal momenta is now dominated by SV, not 
by 0J vac . Since \SV\ » D, the collisions are rare relative 
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to the oscillation period.^ Averaging over an oscillation 
period, and if we begin with one flavor, the average prob- 
ability for the appearance of the other is \ sin 2 (20) where 
is the in-medium mixing angle. If the oscillations are 
interrupted with a rate r co n, we have 



Tfia 



sin 2 (26)L> 
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if D is interpreted as r co n/2. For maximally mixed neu- 
trinos, the in-medium mixing angle is given by 



tan(26) 



SV 



(5) 



Since |u>va C /6V'| <C 1 we have tan(29) w sin(20) so that 
we recover Eq. (|TJ). A large refractive energy difference 
between the flavors suppresses the conversion rate, an 
effect which evidently can be interpreted as a suppression 
of the mixing angle in the medium. 



III. RATE OF DECOHERENCE 

In order to estimate the rate of flavor conversion we are 
thus left with the task of estimating the damping rate, 
D, and the refractive effect SV for the conditions of a 
SN core. The largest possible conversion rate obtains for 
SV = so that we first consider this case. If rfl avor were 
slow in this limit, a further discussion of refractive effects 
would be unnecessary. 
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the diffusion rate r. 
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where we have used Eq. (|^). Evidently it is the low- 
energy neutrino modes with their fast oscillation fre- 
quency which are most effective for flavor conversion. 

There are two conceptually different contributions to 
D. Oscillations are interrupted by scattering processes 
which are sensitive to the neutrino flavor such as + 
e~ — > v e + fJ,~ , a process which has no analogue for 
v T because of the large r mass. In addition, a mixed 
neutrino can scatter to higher energies, for example by 
v+e~ — > e~ +z/, without interrupting the oscillation pro- 
cess. However, after this "up-scattering," the probability 
for processes which do distinguish between the flavors 
increases significantly since the cross sections for all neu- 
trino processes increase with energy. Moreover, at large 
energies the fast beta process + n — > p + fi~ becomes 
kinematically allowed. Effectively, the up-scattering of a 



'Sometimes this situation is referred to as "weak damping." 
However, we follow the convention of Ref. || where "weak 
damping" means D -C i^vac- 



low-energy neutrino causes the oscillations to be inter- 
rupted, even if the up-scattering process itself is flavor- 
diagonal. Therefore, up to a numerical factor, the deco- 
herence rate D of low-energy neutrinos is identical with 
their scattering rate to higher energies. 

The most important energy-changing process for neu- 
trinos in a SN core is the scattering on the nuclear 
medium. Ignoring the subdominant vector-current inter- 
action, the differential cross-section for neutral-current 
scattering of a neutrino of energy E\ to energy Ei is 



da 
dE~2 



3C A G\ 2 S(E 1 - E 2 ) 
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Here, S(u>) is the dynamical structure function for the 
axial-vector current interaction in the long-wavelength 
limit pTJ]. In a dilute medium S(lj) = 2ttS(uj) : lead- 
ing to the usual neutral-current elastic scattering cross 
section. In a nuclear medium, nucleon-nucleon interac- 
tions cause S(u>) to be a broadly smeared-out function 
which obeys the detailed-balancing condition S(— ui) = 
S(u))e~^/ T . Unless spin-spin correlations and degener- 
acy effects are important, the structure function has the 
norm J S(uj)duj /2ir = 1. 

In contrast with elastic neutrino scattering processes, 
the cross section Eq. (fjj) does not vanish for small neu- 
trino energies. With E\ = 0, the neutrino scattering rate 
on nucleons is 



= (0) 
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where rig is the baryon (nucleon) number density. With 
Ca = 1-26/2 for neutral-current processes, the coefficient 
before the integral is 



0.98 x 10 s s~ 



3 x 10 14 g cm 



13 J 30' 



(9) 



where T 30 = T/30 MeV. 

The dimensionless integral vanishes if S(u>) is very nar- 
row, corresponding to the usual case of a vanishing elastic 
neutrino scattering cross section at E\ = 0. The inte- 
gral also vanishes when it is very broad because S(lo) 
decreases at least exponentially for large negative ui. If 
S(ui) is normalized, and if it is a smoothly varying broad 
function, the maximum possible value for the integral ex- 
pression is about 0.25 which obtains when the width of 
S(li) is of order the temperature, probably correspond- 
ing to realistic SN conditions. The rate gets reduced 
by nucleon degeneracy effects and spin-spin anticorrcla- 
tions. Therefore, for low-energy neutrinos the upscatter- 
ing rate is probably not larger than a few times 10 7 s _1 . 
For low-energy neutrinos, this rate is larger than any 
other up-scattering process that we could identify such 
as neutrino-electron or neutrino-neutrino scattering. 

Below the muon production threshold in beta pro- 
cesses, the most important reactions which distinguish 



between and v T appear to be 
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and 
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Vf_ L + M _ — > e~ + v e . We find that for typical conditions 
in a SN core and for neutrino energies around T the rate 
is at most a few times 10 7 s _1 and thus comparable to 
neutrino nucleon scattering. 

Comparing these numbers with Eq. (|2|) reveals that 
D 3> w V ac- Therefore, the strong damping limit indeed 
applies as had been assumed earlier. 

The actual decoherence rate D is smaller than our es- 
timates of the corresponding collision rates. We have 
already mentioned that D is half the collision rate in 
a situation where one flavor scatters while the other is 
sterile B. Altogether, we believe that 



D < 3 



10 7 s" 



(10) 



is a realistic upper limit for typical SN conditions of p 
3 x 10 14 g cm" 3 and T = 30 MeV, implying 



^flavor "^diffusion ^ 20 A.3 E^q 



(ii) 



Averaging this expression over a thermal neutrino distri- 
bution, we note that (E~ 2 ) = 0.38 T~ 2 so that we finally 
conclude that for typical SN conditions 



^flavor) ^diffusion ^ 10 A3 



(12) 



Near the upper end of the mass range suggested by Su- 
perKamiokande we have Sm 2 = 10" 2 eV 2 or A 3 = 10, 
implying that the flavor conversion rate increases by two 
orders of magnitude relative to our estimate. 

In summary, it appears that the flavor conversion be- 
tween za, and v r in a SN core would be much faster than 
the diffusion time scale of about one second if refractive 
effects could be ignored. 



IV. REFRACTIVE EFFECTS 

Turning to neutrino refraction in the SN medium, we 
begin with the usual lowest-order weak potential, 



SV {1) =V2G F (An M + An„ - An T - An„ 



(13) 



where Gf is the Fermi constant while Arij stands for 
the number density of particle j minus the density of 
anti-particles. Since the total fj, and r lepton number 
trapped in a SN core is zero, we have initially SV^ = 0. 
Therefore, the flavor conversion rate would seem to be 
initially fast until a significant weak potential difference 
has built up. 

However, second-order effects can not be neglected. At 
one-loop level, the neutral-current interactions of and 
v r with the nucleons are not identical because of the dif- 
ferent charged-lepton masses in the loop. The second- 
order energy difference in a normal medium was found 
to be Ell 
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6.3 x 10" 4 eV 



6.61 + y;/3 



3xl0 14 g/cm 3 
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(14) 



where hb is the baryon density and Y n the neutron num- 
ber fraction. With 6.3 x 10~ 4 eV = 9.6 x 10 11 s" 1 we 
find \5V^ | > D > w vac and that even initially the flavor 
conversion time scale (the inverse of rfl avor ) far exceeds 
the one for diffusion. This must be the only example 
where radiative corrections to the neutrino refractive in- 
dex are of direct practical relevance! 

We stress, however, that 8V^ is by no means exoti- 
cally large. Its importance in the present context derives 
from the cancellation of SV^ and from the smallness of 
cj vac . For example, the in-medium mixing angle in the 
present case is 

i.e. second-order refractive effects suppress the mixing 
angle by about eight orders of magnitude! The real sur- 
prise is that 07 V ac, despite its smallness, is large enough to 
cause flavor equilibrium if it were not for the refractive 
suppression of the mixing angle. 

Another second-order contribution to 8V arises from 
the low-energy tail of the M /± and Z° resonance in neu- 
trino forward scattering on other neutrinos and charged 
leptons fl2|| . This term is the dominant second-order 
correction in the early universe where the particle- 
antiparticle asymmetries are small, but in a SN core 
Eq. (^) is more important because of the huge baryon 
density. 

Once diffusion processes begin to build up a net /i or 
r lepton number in the SN core, the first-order refrac- 
tive effect becomes important. Its magnitude is under- 
stood if we calculate the refractive energy shift caused 
by neutrinos with a chemical potential fi u <?C T which 
we express as r\ v — ji v /T. To lowest order we have 
An v = T 3 7] v /6 + 0(tjI) so that 



SVP = 0.074 eV T? 



(16) 



Likewise, for muons we find An M = (T 3 /ir 2 ) /(m^/T) ?y M 
where the function f(m tl /T) is 1.077 for m M /T = 3. 
Therefore, muons provide a similar energy shift. Alto- 
gether, if we use rj as some characteristic chemical po- 
tential for the muons and neutrinos, the first-order en- 
ergy shift is given by Eq. ([ll]), leading to an in-medium 
mixing angle of 



tan(29) = 2.2 x 10" 



A a 



^3oTf 



(17) 



Again, the mixing angle is hugely suppressed unless 
rj is finely tuned to zero. Considering that rfl avor = 



4 



sin 2 (20)D the suppression of the conversion rate is 
quadratic in the small in-medium mixing angle. For 
D = O(10 7 s _1 ) even a tiny value for 77 is enough to 
suppress flavor conversion completely. 

As n and r lepton number builds up it may happen 
that the first and second-order contributions to SV cancel 
for suitable conditions, allowing briefly for a fast flavor 
conversion rate. However, the slightest deviation from 
this condition again suppresses the mixing angle and thus 
quenches any further flavor conversion. Any attempt to 
reach flavor equilibrium by neutrino oscillations is ro- 
bustly self-quenching. 

We have mostly studied typical core conditions in 
the SN. In regions below the neutrino sphere the den- 
sity is about three orders of magnitude less than our 
benchmark value of nuclear density, and the tempera- 
ture may be a factor of 5 smaller than our standard 
figure of 30 MeV. For such conditions, the presence of 
muons is strongly suppressed so that the decoherence 
rate D is much smaller than what we have estimated 
for average core conditions. On the other hand, the in- 
medium mixing-ang le Eq. @ is still very small if we 
use T30 = 0.2 and E30 — 0.2 as characteristic values near 
the neutrino sphere. Therefore, our conclusion that fla- 
vor equilibrium cannot be achieved applies to conditions 
throughout the SN core. 



V. CONCLUSIONS 



quences which need to be explored. In the past, most of 
the attention was focussed on the possibility of detect- 
ing evidence for neutrino oscillations in the astrophys- 
ical context. However, as neutrino oscillations become 
more and more experimentally established, the problem 
of unravelling the neutrino mass matrix may become a 
less pressing astrophysical preoccupation than the reverse 
question: given the experimentally measured neutrino 
parameters, is the effect of flavor violation important or 
ignorable in a given environment? 

We think it is intriguing that the core of a SN is pro- 
tected from the consequences of flavor-lepton number vi- 
olation by the phenomenon of neutrino refraction and 
that second-order effects have to be included. While 
the significance of the weak-interaction potential for the 
resonant enhancement of oscillations in the spirit of the 
MSW effect has been widely acknowledged, the suppres- 
sion of oscillations in a SN core is another consequence of 
neutrino refraction with real and important astrophysical 
ramifications. 
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The chemical equilibration between the /i and r flavors 
in a SN core by neutrino oscillations would be fast on the 
neutrino diffusion time scale if the refractive energy shift 
SV between and v T were small. The usual first-order 
contribution indeed vanishes due to the lack of trapped 
/1- and r- lepton number, but second-order contributions 
are large enough to suppress flavor conversion. 

In the course of neutrino transport to the stellar sur- 
face, \i- and r-lepton number will build up due to the 



differences between the v„ 



diffusion con- 



stants ||. It is conceivable that local conditions can be 
reached where the refractive energy shift SV vanishes to 
all orders. While this cancellation would momentarily 
lead to a fast rate of flavor conversion, the resulting re- 
distribution of \x and r lepton number quickly produces a 
SV so large that the conversion rate is suppressed again. 
The equilibrium condition fi u — [i Vt implies a SV so 
large that it can never be reached, i.e. the flavor conver- 
sion process is inevitably self-quenching. 

For the small neutrino mass differences indicated by 
the experiments, the lack of conversion between v e and 
the other flavors is much easier to understand because the 
large amount of trapped electron-lepton number causes 
SV^- 1 ' to be so large that the in-medium mixing angle is 
easily seen to be vastly suppressed (5) . 

The phenomenon of neutrino mixing and neutrino 
oscillations may have a variety of astrophysical conse- 
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